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Abstract 

Let B"' be an (A'i, c?)-fractional Brownian motion with Hurst index {i = 1,2), and let 
B°'^ and B"^ be independent. We prove that, if + > d, then the intersection local times 
of B"^ and i?"^ exist, and have a continuous version. We also establish Holder conditions for 
the intersection local times and determine the Hausdorff and packing dimensions of the sets of 
intersection times and intersection points. 

One of the main motivations of this paper is from the results of Nualart and Ortiz-Latorre 
(J. Theor. Prohab. 20 (2007)), where the existence of the intersection local times of two 
independent (1, (i)-fractional Brownian motions with the same Hurst index was studied by using 
a different method. Our results show that anisotropy brings subtle differences into the analytic 
properties of the intersection local times as well as rich geometric structures into the sets of 
intersection times and intersection points. 

Running head: Regularity of intersection local times of fractional Brownian motions 
2000 AMS Classification numbers: 60G15, 60J55, 60G18, 60F25, 28A80. 
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1 Introduction 

Let Bq = {Bq{u), u G W} be a p-parameter fractional Brownian motion in R with Hurst index 
7 € (0, 1), i.e., a centered, real- valued Gaussian random field with covariance function 

E [Bl{m)Bl{u2)] = \ + \u^f^ - - u^f-^) . (1.1) 

It follows from Eq. (jl.ip that E[(i?Q'(ni) — B'^iu'i)^^ = \ui — and B^ is 7-self-similar with 

stationary increments. 

We associate with Bq a Gaussian random field B"' = {B'''{u),u € W} in M"? by 



B^u) = {Bj{u),...,B^^iu)), ueRP, (1.2) 
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where Bj, . . . , Bq are independent copies of Bq. B^ is called a (p, g)-fractional Brownian motion 
of index 7. 

Fractional Brownian motion has been intensively studied in recent years and, because of its 
interesting properties such as short/long range dependence and self-similarity, has been widely 
applied in many areas such as finance, hydrology and telecommunication engineering. 

Let = {S°i(s), s G M^i} and = {B'^2{t), t G M^^} independent fractional 

Brownian motions in with Hurst indices ai, 02 £ (0, 1), respectively. This paper is concerned 
with the regularity of the intersection local times of B"^ and B"^ , as well as the fractal properties 
of the sets of intersection times and intersection points. Without loss of generality, we further 
assume ai < 02 throughout this paper. For A''! = A''2 = 1 and ai = 02 = ^, the processes are 
classical d-dimensional Brownian motions. The intersection local times of independent Brownian 
motions have been studied by several authors [see Wolpert (1978a), Geman, Horowitz and Rosen 

(1984) ] and is closely related to the self-intersections (or multiple points) of Brownian motion. The 
approach of these papers relies on the fact that the intersection local times of independent Brownian 
motions can be seen as the local times at zero of some Gaussian random field. For the applications 
of the intersection local time theory for Brownian motions, we refer to Wolpert (1978b) and LeGall 

(1985) , among others. 

The self-intersection local times of fractional Brownian motion were studied by Rosen (1987) 
for the planar case, and by Hu and Nualart (2005) for the multidimensional case. Very recently, 
Nualart and Ortiz-Latorre (2007) proved an existence result for the intersection local times of two 
independent d-dimensional fractional Brownian motions with the same Hurst index. 

The aim of this paper is to show that the existence of the intersection local times for two 
independent fractional Brownian motions B'^^ and B°'^ in can be studied by using a Fourier 
analytic method and, moreover, this latter method can be applied to establish the joint continuity 
and sharp Holder conditions for the intersection local times. Besides their own interest, these 
results are useful for studying fractal properties of the set of intersection times as well as the set of 
intersection points. 

Let X = {X{s, t), {s, t) e M^} be an (A^, d)-Gaussian random field, where N = Ni +N2, defined 

by 

X{s,t) = B'^^{s)-B''^{t), sGM^i, tGM^2_ ^1,3) 

We will follow the same idea as Wolpert (1978a) and Geman, Horowitz and Rosen (1984) and 
treat the intersection local times of B"^^ and i?"^ as the local times at of X, with an intension 
to establish sharp Holder conditions. The main ingredients for proving our results are the strong 
local nondeterminism of fractional Brownian motions, occupation density theory [cf. Geman and 
Horowitz (1980)], and newly developed techniques for anisotropic Gaussian random fields [cf. Ay- 
ache, Wu and Xiao (2008) and Xiao (2009)]. 

For later use, we mention that, by the self-similarity and stationarity of the increments of B°^^ 
and -B°2 , the Gaussian random field X defined by ()1.3p has stationary increments and satisfies the 
following operator-scaling property: For every constant c > 0, 

{X{c^ {s, t)), {s, t)ER^}^ {cX{s, t), {s, t) G M^} , (1.4) 

where A = {aij) is an A^ x A^ diagonal matrix such that ajj = 1/ai if 1 < i < A'"! and an = 1/02 if 
A'^i -f 1 < i < A^. In the above, = denotes equality of all finite dimensional distributions and is 
the linear operator on defined by = Yl^=o ^^"'^^i^ • 



2 



This paper is organized as follows. In Section [2l we give several lemmas which will be used to 
prove our main results in the following sections. In Section [3l we study the existence and the joint 
continuity of the intersection local times of two independent d-dimensional fractional Brownian 
motions. We prove that the necessary and sufficient condition for the existence of an intersection 
local times in L^(P x A^^) actually implies the joint continuity. We devote Section S] to the study 
of the exponential integr ability and Holder conditions for the intersection local times. The later 
results imply information about the exact Hausdorff measure of the set of intersection times of B"^^ 
and -B"^ . Finally, in Section [5l we determine the Hausdorff and packing dimensions of the set of 
intersection points of B°^^ and B"'^. 

Throughout this paper, we use (■, •) and | • | to denote the ordinary scalar product and the 
Euclidean norm in R^, respectively, no matter what the value of the integer p is, and we use \p 
to denote the Lebesgue measure in R^. We denote by Op{u,r) a p-dimensional ball centered at u 
with radius r, and Op^^p^(u,r^ := Op^{ui,r) x Op^{u2.,r), where u = {ui,U2) with ui € R^^ and 
U2 € R''^. In Section i, unspecified positive and finite constants will be numbered as c. j, c. 



2 Preliminaries 

In this section, we provide necessary preparations for the proofs of our main results in the later 
sections. 

It follows from Lemma 7.1 of Pitt (1978) that, for any 7 G (0, 1), the real-valued fractional 
Brownian motion Bq = {Bq{u),u E R^} has the following important property of strong local 
nondeterminism: There exists a constant < Cj ^ < 00 such that for all integers n > 1 and all 
u, ui, . . . ,Un e RP, 

Vav(B^{u)\B^{m),...,B^{un)) >c,, min \u - Uk\^\ (2.1) 

where Var(i?Q (ti)|i?Q(ni ),..., i?Q(nn)) denotes the conditional variance of Bq{ u) given i?g('Ux), 
. . . ,i?J(n„), and where uq = 0. The strong local nondeterminism has played important roles in 
studying various sample path properties of fractional Brownian motion. See Xiao (1997, 2006, 
2007) and the references therein for further information. It will be the main technical tool of this 
paper as well. 

We consider the real-valued Gaussian random field Xq = {Xo(s, t), (s, t) G R^} defined by 
Xo{s,t) := Bq^{s) - B^'^it) for s G R^^ and t G R^^. Then the coordinate processes of X defined 
by (jl.3p are independent copies of Xq. 

The following Lemma 12.11 is a consequence of the property of strong local nondeterminism of 
fractional Brownian motion, and will be useful in our approach. 

Lemma 2.1 There exists a constant < c^^ < 00 such that for all integers n > 1 and all 
{v,w), (si,ti), . . . , {sn,tn) G R^, wc havc 

Var(Xoiv,w)\Xo{si,ti),. . . ,Xo{sn,tn)) >c,J min \v - Sk]"^"' + min \w - tkl'^'^A , (2.2) 

\ V 0<k<n 0<k<n J 

where sq = to = 0. 
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l<i<n 



Proof By definition we can write 

Yai(Xo{v,w)\Xo{si,ti),...,Xo{sn,tn))= inf E 

\ / aiGK, 1 

Since B^^ and Bq'^ are independent, we have 
Vav(^Xo{v,w)\Xo{si,ti), . . . ,Xo(s„,t„)^ 

E 



Xo{v, w) - '^ aiXo{si,ti) j 
i=i ' 



(2.3) 



inf 

aiGK, l<i<n 



> inf E 

aiGM, \<i<n 



n . 

B^^{v)-Y^a,B^^{sM 
1=1 

n 

B^^{v)-Y,a^B^^{s.) 



+ E 



1=1 



+ inf E 

feiGK, l<i<n 



n 

Br{w)-Y,hBr{U 



1=1 



(2.4) 



= YaT[B^^{v)\B^^{si), . . .,B^^{sn)) + Xavi^B^^w)\B^^h), . . .^B^^t^) 
Hence (IHSD follows from ((231) and (lO. 



Combining Lemma 12.11 with the following well-known fact, which will be used repeatedly 
throughout the paper, that 



detCov(Zi, ...,Zn) = Var(Zi) Var(Zfc|Zi, . . . , 



(2.5) 



k=2 



for any Gaussian random vector (Zi, . . . , Z„), we have that, for any (si, ti), . . . , {sn,tn) £ , 



detCov (Xo(si, ti), . . . , Xo{Sn, tn)) 



an 



Var(i3o"^(.,)|i3-^(.i), . . . ,i?o"^(s,-i)) + Var(i3-^(t,)|i3o"^(ti), . . . , B^^t,^^)) 



(2.6) 



To prove the existence and continuity of the intersection local times of B"^ and , we will 
make use of the following lemmas. Lemma 12.21 is similar to Lemma 8.6 in Xiao (2009) whose proof 
is elementary. Lemma 12. 3 1 and Lemma 12.41 are extensions of Lemma 2.3 in Xiao (1997) and will be 
useful for dealing with anisotropy of the Gaussian random field Xq . Lemma 12.51 due to Cuzick and 
DuPreez (1982), is a technical lemma. 



Lemma 2.2 Let (3, 7 and p be positive constants, then for all A G (0, 1) 



1 ^p-i 



(A + r'r)' 



dr >i < 



log (1 + yl-l/T) 



if I3j > p, 
if j3j = p, 
if l3j < p. 



(2.7) 



In the above, f{A) x g{A) means that the ratio f{A)/g{A) is bounded from below and above by 
positive constants that do not depend on A £ (0, 1). 
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Proof This can be verified directly and we omit the details. □ 

Lemma 2.3 Let (3, 7 and p be positive constants such that 7/3 > p. 

(i). If > p, then there exists a constant 3 > whose value depends on 7, /3 and p only such 
that for all A G (0, 1), r > 0, G M^, all integers n > 1 and all distinct ui, . . . ,Un € Op{u* ,r) 
we have 

du P-3 

g<C2^nA'' (2.8) 

Op(u*,r) [A + min{|ti — UjY^ , j = 1, . . . ,n}f 

(a). If = p, then for any k £ (0,1) there exists a constant c^ ^ > whose value depends on 
7, (3, K and p only such that for all A (0,1), r > 0, u* £ W^, all integers n > 1 and all 
distinct ui, . . . ,Un G Op{u* ,r) we have 



du 

<c^,n los 



1/3 - --2,4 



lOp(w,r) [A + min{|'u - Uj\^, j = 1, . . . ,n}] 

Proof The idea of proof is similar to that of Lemma 2.3 in Xiao (1997). Let 
Lj = < n G Op{u*,r) : \u — Ui\ = min{[u — Uj\, j = 1, - ■ ■ ,n} 



e + 



(2.9) 



Then 

n n 

Op{u*,r) = \Jri and Xp{Op{u* ,r)) =^ Xp{Ti). (2.10) 

1=1 i=l 

For every u € Lj, we write u = Ui + p9, where < p < Pi{9) and 9 G the unit sphere in W. 

Then 



r rpr{0) 
\{Ti) = Cp / u{d9) / ff-^dp 

= ^ / prn'^^ido), 

where v is the normalized surface area in Sp-i and Cp is a positive constant depending on p only. 

Denote the integral in (j2.8p and (j2.9p by Ii. We first consider the case of 7/? > p. By (j2.10p . a 
change of variables and Lemma 12.21 we can write Ii as 



du 



1=1 

n 



■[Jr, [A + mm{\u - Ujl'y , j = l,...,n}f 

Sp-1 ' 'Jo {A + p-yy 



V- ^31 , . fA-'Z-'pm pP-l (2.12) 

Va^7^^ / u(d9) / — -^dp ^ ' 



i=l •'Jp-i 



n „ 

<c,,,Y,A'^-n u{d9) 

A 1 Srt^l 



C2 3 n A 



i=i -J^p-i 
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This proves inequality (j2.8p . 

Now we assume 7/3 = p. As above, we use (|2.10p and a change of variables to get 



< — E / log 

i=l 



(l + p7)/3 



dp 



e+(A-'/"'p,ie) 



(2.13) 



In the above, we have used the fact that if 7/3 = p and k E (0, 1), then for all x > 

d/9 < -log(e + x''). 



Since the function tlJi{x) = log(e + x"/^) is concave on (0,oo), we apply (12.111) and Jensen's 
inequality twice to derive 



n 

/i<c,,3n^-V'i(vl-P/%(ri) 

i=l 



<c^^n log 



e + 



This finishes the proof of (|2.9p . 



(2.14) 



□ 



Lemma 2.4 Let P > be a constant and let p > 1 be an integer such that (3 < p. Then the 
following statements hold: 



(i). For all r > 0, u* £ M^, all integers n > 1, and all distinct ui, . . . ,Un £ Op{u*,r), we have 

du 



I 

J Op{u* ,r 



•j min{[u — Uj\f^, j = 1, . . . ,n} 



< c^f^np r 



(2.15) 



where c^^ > is a constant whose value depends on (3 and p only. 



(a). For all constants r > and K > 0, all u* G M^, integers n > 1, and all distinct ui, . . . G 
Op{u*,r), we have 



I 

JOJu* 



log 



p{u*,r) 



< c^^r^ lo, 



e + K(^mm{\u — Uj\, j = 1, . . . , n}) 

r \-0 



-13 



du 



e + K 



n 



i/p 



(2.16) 



where c^^ > is a constant whose value depends on (3 and p only. 
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Proof Part (i) is a special case of Lemma 2.3 in Xiao (1997). Hence, it only remains to prove 
Part (ii). Denote the integral in (j2.16p by I2. As in the proof of Lemma 12.31 we have 



I2 / log e + K { mm{\u - Uj\, j = 1, . . . , n}) 



du 



(2.17) 



= Y.Cp v{d9) / log [e + K p-^) dp 

j—]^ Sp-i Jo 

= Y.^p p'^^y ^(^^) / f^'^ + Kpm-Cp-") dp 

1—1 J Sp-\ Jo 
n 

<c,,s Y. / P^{0)nog{e + Kpm~^)v{de). 

i=l J Sp-1 

In deriving the last inequality, we have use the fact that log(e + xy) < log(e + x) + log(e + y) for 
all X, y > 0. Since P < p, we can verify that the function ^2ix) = x log{e + K x-^/P) IS concave on 
(0,00). By using Jensen's inequality twice, we obtain 



h < s X] ^2 ( / 

i=l ^ •S'p-i 

/ 1 " \ /'rP 

^ i=l ^ 



(2.18) 



This finishes the proof of Lemma 



□ 



Lemma 2.5 Let Zi, . . . , Z„ he the mean zero Gaussian random variables which are linearly inde- 
pendent and assume that 

coo 

g{v)e~^'" dv < 00 



for all e > 0. Then 



g{vi)exp --Var( VwjZj 

(27r)("-i)/2 
(detCov(Zi,...,Z„))V2 



dvi ■ ■ ■ dVn 



^)e-^y'dv, 



(2.19) 



where a\ = Var(Zi|Z2, . . . , Z„) is the conditional variance of Z\ given Z2, . . . , Z^. 

3 Intersection local times and their joint continuity 

In this section, we briefly recall the definition of local time as occupation density [cf. Geman and 
Horowitz (1980)] and then study the existence and joint continuity of the intersection local times 
of and 
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Let Y{t) be a [random] Borel vector field on W with values in M"^. For any Borel set E (^MP, 
the occupation measure of y on S is defined as the following measure on W: 

If fi^ is absolutely continuous with respect to the Lebesgue measure Xg, we say that Y{t) has 
local time on E, and define its local time, L{»,E), as the Radon-Nikodym derivative of fi^ with 
respect to Xg, i.e., 

dXg 

In the above, x is the so-called space variable, and E is the time variable. Note that if Y has local 
times on E then for every Borel set F C E, L(x, F) also exists. 

It follows from Theorem 6.4 in Geman and Horowitz (1980) that the local time has a measurable 
modification that satisfies the following occupation density formula: For every Borel set E CMP, 
and for every measurable function / : R'' ^ M-(-, 



f{Y{t))dt= / f{x)L{x,E)dx. (3.1) 

Suppose we fix a rectangle E = [a, a + h] Q W, where a E and h € MP^. If we can choose a 
version of the local time, still denoted by L{x, [a, a + t]), such that it is a continuous function of 
{x,t) € W X [0,h], Y is said to have a jointly continuous local time on E. When a local time is 
jointly continuous, L{x, •) can be extended to be a finite Borel measure supported on the level set 

Y£\x) = {teE:Y{t)=x}; (3.2) 

see Theorem 8.6.1 in Adler (1981) for details. This makes local times, besides of interest on their 
own right, a useful tool in studying fractal properties of Y. 

It follows from (25.5) and (25.7) in Geman and Horowitz (1980) that, for ah x,yeRi,EQ W 
a closed interval and all integers n > 1, 



E 



Lix,Er] = (2vr)-"^ / / exp f -if;(n„x)') 

X Eexp ^ {uj,Y{tj))^ du dt 



(3.3) 



and, for all even integers n > 2, 

E 



(L(x, E) - L{y, S))"l =(27r)-"'' / / TT fe-^<"^'^> - e'^^"^'? 

X Eexp ii'^{uj,Y{tj)) ]dlL dt 

^ 7=1 ^ 



(3.4) 



where u = (ui, . . . , Un), t = (ti, . . . , tn), and each uj G M'^, tj € E. In the coordinate notation we 
then write Uj = {uj^i, . . . , uj^g). 

The main results of this section are the following Theorem 13. II and Theorem 13. 31 for the existence 
and the joint continuity of the intersection local times of two independent fractional Brownian 



motions in M*^. 
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Theorem 3.1 Let B"^ = {^"^(s), s E M^i} and = {B°"'{t), t G R^^} be two independent 
fractional Brownian motions with values in R'* and Hurst indices ai and a2, respectively. Let 
X = {X{s,t), {s,t) S M^} be the {N,d)- Gaussian random field defined by liL^J]) . Then, for any 
given constant R > 0, X has a local time L{x, Oni,N2{^i B)^ G L^(Px A^) if and only ^ + ^ > d. 
Furthermore, if it exists, the local time of X admits the following L"^ -representation 

L{x,ON,,NMR)) = {'^^r'' I e-*<^'^> / e^<^'^"^W-^"^W>dsdtdy, (3.5) 



Vi,jV2 > 



and the local time L can be chosen as a kernel L(-, •) on X B{Oni,N2{0, R)) ■ In particu- 
lar, if — -\- — > d, then B"^ and B°^'^ have an intersection local time which can be defined as 
L"i'"2 {ONuNTiO, R)) ■=L{0, On.MO, R)) ■ 

Some remarks about Theorem 13. II are in order. 

Remark 3.2 (i) When Ni = N2 = 1, ai = a2 = H and Hd < 2, the existence of the intersection 
local time was proved by Nualart and Oritz-Latorre (2007) as the L'^-limit of 

Ie{B^,B")= / pe{B^is)-B"{t))dsdt, as e^O, (3.6) 
Jo Jo 

where Pe{x) = (27re)'^/^ exp(— |xp/(2e)). They also proved that if Hd > 2, then 
limE\UB^,B")] = 00 and lim Var [4(5^, 5^)1 = 00. 

elO ^ ^ eiO ^ ^ 

In the above, B^ = {B^{t),t > 0} and B^ = {B^{t),t > 0} are two independent fractional 
Brownian motions with values in and index H € (0, 1). Similar method can be applied to show 
that the intersection local time L"!'"^ ((^^^^^^(o, i?)) in Theorem 13.11 can be chosen as the L^-limit 
of the following approximating functionals 

= / pe{B'''is)- B''%t))dsdt, as e^O. (3.7) 

JOn^.N2{0,R) 

Moreover, we are able to show that, if — + — = d, then 

' ' ai 02 ' 

E[/,(5°S5°2)] ^c(ai,a2,iVi,A^2)ln0^ , as e ^ 0, (3.8) 

where c{ai,a2, Ni, N2) > is a constant depending on ai, 02 and A^i, N2 only. This raises an 
interesting question whether Ig can be renormalized to converge to a non-trivial limiting process. 
This and other related questions will be dealt with elsewhere since they require different methods. 

(ii) It follows from the operator-scaling property (|1.4p of X and (j3.5p that the intersection local 
time ((9^^^^2(0, i?)) has the following scaling property: For any constant c > 0, 

{L"i'"2(c^0;vi,^2(0,i?)), R>0} = {c^+^"^L-i'"2(O^„^2(0,i?)), R>0]. (3.9) 
Here A is the N x N diagonal matrix as in (jl.4p . 
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(iii) We say that the sample functions of B'^^ and B°'^ intersect if there exist s € and 
t € such that B°'^{s) = It is also of interest to study the geometric properties of the 

set of intersection times 

M2 = {{s,t) G : = B'^'it)} 

and the set of intersection points 

D2 = {x eR'^ ■.x = B°'^{s) = for some {s,t) G R^}, 

because they are often random fractals. The existence of the intersection local time and its prop- 
erties are closely related to the existence of intersections of the sample functions of B°'^ and 
and the geometric properties of M2 and D2- Similar to Theorem 7.1 in Xiao (2009), we can prove 
that if + > then M2 7^ with positive probability. On the other hand, Theorem 3.2 in 
Xiao (1999) proved that if ^ + ^ < d then M2 = almost surely. In Section H we will give more 
information on the Hausdorff and packing dimensions of M2 , as well as a lower bound for the exact 
Hausdorff measure of M2. The Hausdorff and packing dimensions of D2 are determined in Section 

Proof of Theorem 13.11 Note that the Fourier transform of the occupation measure u„ 
of X is 

It follows from the Plancherel Theorem that X has a local time -L(x, O7Vi,Af2(0) -^)) £ L'^{^ x A^) 
with a representation (13. 5p if and only if 

J := / dsdtdvdw Eexp (i{y, X{s,t) - X{v,w))) dy < 00. (3.10) 

See Theorem 21.9 of Geman and Horowitz (1980). Hence, it suffices to prove that Eq. ()3.10p holds 
if and only if ^ + g > d. For this purpose, we use the independence of the coordinate processes 
of X, (fO|) andl^J^i^to deduce that 



J 



dsdtdvdw 



E{Xo{s,t)-Xo{v,w)Y 



dsdtdvdw 



(3.11) 



By using spherical variable substitutions and Lemma 12.21 it is elementary to verify that the last 
integral in Eq. (|3.11|) is finite if and only if + > d. 

When the later holds, one can apply Theorem 6.3 in Geman and Horowitz (1980) to choose 
a version of the local time of X, still denoted by L, such that it is a kernel in the following 
sense: For every x € M'^, L{x,-) is a finite measure on B{Oni,N2{^i I^)) and, for every Borel set 
E G i3(OAr^^7V2(0, -R)), X I— > L{x,E) is a measurable function. This proves the main conclusion of 
Theorem 13.11 Finally, by taking a; = we prove the last conclusion of Theorem 13. 1[ □ 
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c 

B"'^ have almost surely a continuous intersection local time on ^ 



Theorem 3.3 Let and defined as that in TheoremlSJl V ^ + ^ > d, then and 



As in the proof of Theorem 13.11 we will prove a stronger result that X has almost surely a 
jointly continuous local time on M^^^^^^ proof is based on the following Lemma 13.41 and 

Lemma 13.61 They will also play an essential role in Section H] for establishing Holder conditions for 
the intersection local times. 

Under the condition ^ + ^ > d, define 

_ / 1 if ^ > d, 

( N-aid if r = 1, 

N2 + iVi - a2d if r = 2. ^^^^> 

[Recall that we assumed ai < 02 throughout the paper, and = A^i + N2.] We will also make use 
of the following notation: 

if T = 1, 



and 




(3.14) 

Note that, if Ni = aid, then Pr = N2 and rjr = 1. To emphasize the importance of /3r and rjr, we 
point out that /3r is the Hausdorff dimension of the set M2 of intersection times and rjr is useful 
for determining the exact Hausdorff measure of M2. See Section [4] for more information. 

Lemma 3.4 Suppose the assumptions of Theorem \ 3.3\ hold. Then, there exist positive and finite 
constants e £ (0, 1/e) and Cg^, which depend on ai,a2, Ni, N2 and d only, such that for all 
r G (0, e), D := 0]\fi,N2{u,r), where u = (ni,M2) ^ , all x and all integers n > 1, we have 



¥.[L{x,DY] < { 



( c^^ (n!)"! r-^i ^ if^,>d, 



!r"-^n-.,log(e+^^5^) ^f^ = d, (3-15) 



3,1 ^ ' •' CKl Qi ' a2 



In the above, = max{y, 0} for every y € 



Remark 3.5 From (I3.13P and (j3.14p . it can be verified that 

^^^^<Vr<N-Pr. (3.16) 

We observe that the power of n\ in (j3.15p becomes (A^ — (3r)/N when X is an isotropic Gaussian 
field as in Xiao (1997) and is N — Pt when X is anisotropic in every coordinate (with the same 
scaling or Holder index) as in Ayache, Wu and Xiao (2008). These seem to be the extreme cases. 
In the present paper, if we assume A"! / A^2 and ai 7^ 02, then strict inequalities in ()3.16p may 
hold and if, in addition, A^i = aid, then extra logarithmic factors appear in the estimate (j3.15p . 
Lemma 13.41 suggests that the local time L{x, •) may satisfy a law of the iterated logarithm which 
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is different from those for the local times of an {N, (i)-fractional Brownian motion or an (N, d)- 
fractional Brownian sheet with index (a, . . . , a); see (|4.10p . (|4.12p and (j4.13p below. This leads us 
to expect that the exact Hausdorff measure function for M2 may be different from those for the 
level sets of fractional Brownian motion and fractional Brownian sheets, respectively. It would be 
interesting to investigate these problems. 

Proof of Lemma 13.41 Even though the proof of Lemma 13.41 follows the same spirit of the proofs 
of Lemma 2.5 in Xiao (1997) and Lemma 3.7 in Ayache, Wu and Xiao (2008), there are some subtle 
differences [see the remark above]. Hence we give a complete proof. In particular, we provide a 
direct way to estimate the last integral in (|3.17p below. We believe that this method will be useful 
elsewhere. 

It follows from ()3.3p and the fact that Xi, . . . ,Xfi are independent copies of Xq that, for all 
integers n > 1, 



E[L{x,Dr] < (27r)-"'^ / IT I / ^^^P f " i^ar f V ^ Xo(s,-, t,)) (Mk } dt 

Jd" Ur" I \~[ J 

= (2^)-nd/2 /■ rdetCov(Xo(si,ti),...,Xo(s„,t„)) 



d 

' dt, 



(3.17) 



where = (ui^k, • • • > ^^n,fc) ^ I^") t = (si, ti, . . . , s„, t„) and the equality follows from the fact that 
for any positive definite n x n matrix F, 

/ ME>i;^exp(-i.'r.)d.= L (3.18) 



In order to prove Eq. (I3.15p . we consider the three cases separately: ^>d, ^<d<^ + ^ 



and ^=d. 

ai 

In the case that ^ > d, thanks to Eq. (|2.6p . we have 



E[L{x,Dr]<c-^J fl 



■ ds dt 



mm 



{\Sj-S^\"^^, 0<i<j-l} 



S,2 



" 1 \ _ 

, , TT w d^ dt (3.19) 

o%^{u2,r) \Jo-^{u,,r) {mm{\sj - Si\'^\ 0<i<j- l}r ' 



n nN2 



l,{u.,r)% min{|s,- - s,|°i^ < i < i - 1} 



where sq := 0, s = (si, . . . , s„) and t = (ti, . . . , t„). 

Since A''! > aid, we integrate the last integral in Eq. ()3.19p in the order dsn-, ■ ■ ■ , dsi and apply 
Part (i) of Lemma 12.41 iteratively. This yields 

E[L(x,L>)"] < c^^(n!)^ ^n{iVi-ai<i) ^ ^nN2 ^ c^^^(n!)''i r"'^!, (3.20) 

which proves Eq. (j3.15p for the case ^ > d [i.e., r = 1]. 
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In the second and third cases [i.e., < d < + we use (|3.17p and ()2.6p to obtain 

1 



E[L{x,Dr]<clf n 



min \sj — Sfcl"^^ + \tj — ifc|"^) 



(3.21) 



0<A:<j-l 



o<fc<j-i 



To estimate the last integral in (j3.2ip . we will integrate in the order of dsn, dtn, ■ ■ ■ , dsi, dti. In the 
case of < d < + we apply Part (i) of Lemma [231 with A = mino<fc<n-i \tn — to 
derive 



dsr 



ON,{ui,r) ( min \sn — Skl"^ + min \tn — tkl^'^Y 

^0<fc<n-l 0<fc<n-l ^ 



< 



[ min \tn-tkM'^ 

'-0<fc<n-l 



Since a2(d - ^) < iV2, it follows from (13:22]) and Part (i) of Lemma [23] that 



dSfidtyi 



D ( mm \Sn 

0<fc<n-l 



< C2 , n 



min |t„-ifch2)' 

0<A;<n-l ' 
dtn 



^0<fc<n-l ^ 



(3.22) 



(3.23) 



a2 



Repeating the above procedure yields (j3.15p for the case of < d < + 

Finally, we consider the case of ^ = d. Let k G (0, 1) be a constant such that ^02/01 < -/V2. 
Applying Part (ii) of Lemma 12.31 with A = mino<fc<n-i \tn — ifel"^) we have 



dsr. 



ON-,{ui,r) { min |sn — Sfcl"! + min Un — ifcl"^)' 

^0<A:<n-l 0<fc<n-l 



(3.24) 



< nlog 



e + ( ( min \tn - ^ 



It follows from (|3.24p and Part (ii) of Lemma flM (with (3 = «;a2/ai and K = {rn 1/^1)'^) that 

dSrjridtn, 



d( min |s„-Sfc|°i+ min \tn - tkl^'^Y 

^0<k<n-l 0<k<n-l ^ 



< C,^ n 



log 



Ojv2(M2,r') 



< C3 7 n log 
= C3 7?7, log 



e + 



e + 



e+( min |t„-tfc|) 

-Ka2/ai 



-K02/O!l 



n 



l/Afi 



1/Af2 



n 

77, iVi \ k/qi 



(3.25) 



< c, „ n log ( e + 



f.a2-ai 
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[Recall that = max{y, 0}.] 

By iterating the procedure and integrating dsn-i, dtn-i, . . . , dsi,dti, we obtain that 



E [L{x, D^] < (n!)''2 r"^^ "Q i^g j ^ + 
This finishes the proof of the moment estimate (j3.15p . 



(3.26) 



□ 



The following lemma estimates the higher moments of the increments of the local times of 
X. Combined with Kolmogorov's continuity theorem, it immediately implies the existence of a 
continuous version of x i— > L{x, D). 

Lemma 3.6 Suppose the assumptions of Theorem \3.3\ hold. Then, there exist positive constants 
Cgg and Ki, depending on e,ai,a2, Ni, N2 and d only, such that, for any r > 0, D := 0^^^i\:[^{u,r) 
for u = (ui,U2) € M^, all x, y (zW^ with \x — y\ < 1, all even integers n > 1, and all 7 € (0, 1) 
small enough, we have 



E [{L{x,D) - £»))"] < c"^^ (n!)''-+^i^ \x - ^n(/3.-Ki7), 



(3.27) 



Proof Let 7 £ (0, 1) be a small constant whose value will be determined later. Note that by the 
elementary inequalities 

_ ^[ < 2i-T|n|T foralluGM (3.28) 
and \u + v\'^ < \u\'^ + \vp , we see that for all ui, . . . , Un, x, y G M.'^, 



n 



-i{uj,x) _ -i{uj,y) 



< 2(1-^)" \x - Yl 



where the summation ^ ' is taken over all the sequences {ki, . . . ,kn) G {!,..., d}". 
It follows from (j3.4p and (j3.29p that for every even integer n > 2, 



E 



(Lix, D) - L{y, D)Y < (27r)-"^2(i-^)" |x - y^ 



n 

771=1 



Eexp I — i'Y^ -^i^j^^j)) ] dudt 



dt 



n 

m=l 



\um,kmr'^ exp 



1 / " 

-Var( ^iuj, X{sj,tj)) 



l/n 



du 



(3.29) 



(3.30) 



where the last inequality follows from the generalized Holder inequality. 

Now we fix a vector k = {ki,k2, ■ . ■ , kn) €{!,..., d}" and n points {si,ti), . . . , (s^, tn) G -f^\{0} 
such that si, fi, . . . , s„, t„ are all distinct [the set of such points has full nA^-dimensional Lebesgue 
measure]. Let M. = M.{k,t,^) be defined by 



\um,kmr^ exp 



^Var(^f;(n,-, X(s,,t,))) 



l/n 



du 



(3.31) 
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Note that Xi (1 < £ < d) are independent copies of Xq. By the strong local nondeterminism of 
fractional Brownian motions Bq^ and Bq^ and Eq. (j2.6p . the random variables Xi{sj,tj) (1 < ^ < 
d,l < j < n) are linearly independent. Hence Lemma 12.51 gives 



\uT r exp 



1 / " 

-Varl ^{uj,X{sj,tj)) 



(27r)("^-i)/2 



[detCov(Xo(si, ti), . . . , Xo(s„, i„))] 



d/2 



— I e 2 dv 

On 



(3.32) 



< 



n7 1 



[detCov(Xo(si, ti), . . . , Xo(s„, 

where cr^ is the conditional variance of Xfc^(sm, im) given Xi{sj,tj) {i ^ km ot i = km but j ^ m), 
and the last inequality follows from Stirling's formula. 
Combining (jS.Sip and (|3.32p we obtain 



M < 



[detCov(Xo(si,ti), . . . ,Xo(s„,t„))] 



d/2 



n 



(3.33) 



m=l 



The second product in (j3.33p is a "perturbation" factor and will be shown to be small when 
integrated. For this purpose, we use again the independence of the coordinate processes of X, (j2.2p 
and (12.11) to derive 



(3.34) 



C7^ = \ax(^Xk^{Sm:tm) Xk^{sj,tj), j ^ m) 

> "^^^{BZ^lBZis,), j / m) +Y^r[BZ{tm)\BZit,), J + rn 
- 4i2 ~ • i / +min{|tm - . j + rn 

As in the proof of Eq. (I3.15p . we will prove Eq. (I3.27P by cases. 
If ^ > d, then we take 7 G (O, ^(^ - d)) so that 

ai(d + 27) < iVi. 



(3.35) 



For any n points (si, ti), . . . , (s^, t^) € D\{0}, we define a permutation vr^ of {1,2,..., n} such 
that 

I^TT^i)! =min{|si|, i = l,...,n}, 

\s-K,(j) - s,r4j-i)l = min{|si - s,r4i-i)|, i G {1, • • . , n}\{7r3(l), . . . ,7rs(i - 1)}} . 
Then, by ()3.34p . we have 

n 



m=l m= 



^ s:r. n 
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< c: 



■-|^ l''7rs(m) •'7rs(m-l) 
1 



|2ai7 



(3.36) 



[detCov(i3o"^(.i),...,i?o"Hsn))]^' 
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It follows from (f333|) . I^M) and (|336]) that 



3,11 v^'i ' TT ^ 

ii rr?! 



M < 

[detCov(5o"i(si),...,i?o"^(5„))] m=l 



d/2 11 ^7 



[detCov(BQi(si), . . .,BQ\sn))y 

n 

- 3,14 V 11 min{|sj -Si|°i(^+27),0 < i < j - 1}' 



Therefore, by (13.350 and Lemma 12.41 we have 



[ X(fc,t,7)(it < c" , (n!)^ / ft ^-7J7T-^ T dt 

Jon ^ ' - ^ ^ Jd^ fJi min{|s, - Si|"i('i+27),0 < i < j - 1} 



< c 



•^3,15 



(j^!)^i + (^+l)T^n(/3i-2ai7)^ 

We combine ([QO]) and (lOHIl to obtain 



E 



(L(x, i?) - L{y, D)y] < cl^^ (n!)''^ + (^+')^ \x - y^^ r'^i/^i-^^n) , (3.39) 



By choosing the constant ki > max{^^ + l,2ai}, we prove Eq. (13.270 for the case — > d [i.e., 
r = l]. 



.... .coc ^^^^^^^^^ 

Wu and Xiao (2008), we choose 



Now we prove Eq. p.27p for the case of-^ < d < tj^ + t^- Inspired by Lemma 3.4 in Ayache, 



7 e 



4 V ai 02 / / 



and set 



Clearly, we have 



and 



1 ( ,Ni N2 ] , , 
-min<l, ai 1 d),aij> (3.40) 

2 [ ^ai 02 J 

i = ^, i = l-l (3.41) 

pi>l, P2>1, — + — = 1 (3.42) 

Pi P2 



^ = N^-6<N,, 
Pi 

a2d / , ^1 , 

P2 V Oil a\, 



(3.43) 
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where the last inequahty follows from the fact that 6 < ^ (-^ + — . By a simple computation, 
^ + ^ = N, + a,d-^N, + (^-l)5 = N-P, + (^-l)5 (3.44) 

and 



Nipi N2P2~ N2 aiN2 \a1N2 Ni) "'^ VaiA^2 Ni) ' ^' ' 

Furthermore, from the way we define 7, 5 and p2-, we know 

— + 2a27<A^2- (3.46) 

For any n points (si, ti), . . . , {sn-, tn) G -C)\{0}, we define a permutation tt^ of {1,2,..., n} such 
that 

1*714(1)1 = min{|ti|, i = 1,. . . ,n}, 

I^TTtO') - *7rt{i-i)l = min{|ti - i G {1,. . . , n}\{7r4 (1), . . . ,7rt(j - 1)}} . 

Then, by ()3.34p . we have 



n ^ n 



1 



- ^'^^ [detCov(i?o"^(ti),...,,i?o"^(t„))]^- 

Recall from (fTel) that 

detCov (Xo(si,ii), . . . , ,Xo(s„,t„)) 

> detCov {B^' (si), . . . , , So"' (^n)) + detCov {B^^ (h), . . . , , Sq"' (*n)) • 

Hence, 

_i 

detCov (Xo(si,ti ),..., Xo(s„,t„)) ' < detCov {B^'{si),...,B^'isn)) 



2pi 



detCov iB^''ih),...,B^^tn)) 



2P2 



It follows from (1333]) . dOT]) and (|3:i9]l that 



[detCov(Xo(si, ti), . . . , Xo(«„,t„))]"" "1 



^ 3,17 



d/2 



[detCov (5°^ (si), . . . , i?°^(^n))] ^ [detCov (So"^(ti), . . . , , 5o"^(tn)) ] ' 



(3.47) 



(3.48) 



(3.49) 



(3.50) 
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Combining (l330l) . SSMh . S^Mh and Lemma [131 we obtain 



M{k,t,^)dt<c^^^{n\y [ ^ 

' •^oj^iK-'^) [detCov (5^1 (si ),..., 



/ 



3,18 ^ ' 



(3.51) 



In the above, the constant ki > is chosen appropriately by taking into account p.44p . p.45p and 
(I3.39p . The value of ki depends on qi, 02, ^1,-^2 and d only. 
We combine p.30p and (j3.5ip to obtain 



E 



(L(X, LI) < (n!)''2+'^iT|3._y|n7^n(/32-M7). 



(3.52) 



This proves (KTlh for the case of^<ci<-^ + ^ [i.e., r = 2]. The proof of Lemma EJ] is 
complete. □ 

Now we are ready to prove Theorem [37 



Proof of Theorem 13.31 The proof of the joint continuity of the local time of X is similar to that 
of Theorem 3.1 in Ayache, Wu and Xiao (2008) [see also the proof of Theorem 8.2 in Xiao (2009)]. 
Hence we only give a sketch of it. 

It suffices to show that for any fixed u = (ui, U2) G and R > 0, the local time L{x, (s, t)) := 
L{x, [ui, ui+s] X [u2,U2+t]) has a version which is continuous in (x, s, t) G M'^x [0, R]'^ almost surely. 
For simplicity of notation, we assume u = (0,0). Observe that for all x,y £ W^, (s,t), {v,w) G 
[0, i?]^ and all even integers n > 1, we have 



E 



[L{x, [0, s] X [0, t]) - L{y, [0, v] x [0, w])Y 



< 2 



n-l 



+ E 



E (L(x, [0, s] X [0, t]) - L{x, [0, v] x [0, w]))"' 
L{x, [0, v] X [0, w]) - L{y, [0, v] x [0, H))" 



(3.53) 



Since L{x,-) is a finite Borel measure, the difference L(x, [0, s] x [0, t]) — L{x,[0, v] x [0, w]) 
can be bounded by a sum of finite number of terms of the form L{x,Dj), where each Dj is a 
closed subset of [0, R]^ of the form ONii-,r) x 0]\i^{-,r) with the radius r < ^|(s,t) — {v,w)\ := 
^y^|s — + |i — w^P- We can use (j3.15p to bound the first term in (j3.53p . On the other hand, the 
second term in (j3.53p can be dealt with by using (|3.27p . Consequently, there exist some constants 
7 G (0, 1) and uq such that for all x, y G M'^, (s, t), {v, w) G Oni,N2{0j R) (^'^^ and all even integers 
n > no. 



E 



(L(x,[0, s] X [0, t])-L{y,[0, v] x [0, w]))' 
<c^,^{\x-y\ + \{s,t)-{v,w)\y\ 



(3.54) 
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It follows from (j3.54p and the multiparameter version of Kolmogorov's continuity theorem [cf. 
Khoshnevisan (2002)] that there exists a modification of the local times of X, still denoted by 
L{x, [s,t)), such that it is continuous for x E W^, (s,t) G [0, R]^ . This finishes the proof of 
Theorem I3.3i □ 



4 Exponential integrability and Holder conditions for the inter- 
section local times 

In this section, we investigate the exponential integrability and asymptotic behavior of the local 
time L{x, •) of X. As applications of the later result, we obtain a lower bound for the exact 
Hausdorff measure of the set M2 of the intersection times of B"^ and B°^^ . 

The following two technical lemmas will play essential roles in our derivation. 

Lemma 4.1 Under the conditions of Theorem \3.3l there exist positive and finite constants e G 
(0, 1/e), J and c^^, depending on ai, 02, Ni, N2, and d only, such that the following hold: 

anc 



(i) For all {01,02) € and D = Oni,N2 {{0-1^(^2),''') with radius r G (0,e), x G M'^ and all 
integers n > 1, 



E 



L(x + X(ai, 02), D) 



< < 



^n^^l^nN2 n^=ilog (e + 



if^^d, 

if^=d. 



(4.1) 



(ii) For all (01,02) G and D = Oi\]^^N2{{ai,a2),r) with radius r > 0, a;, y G M"^ with 
\x — y\ < 1, all even integers n > 1 and all 7 G (0, 1) small, 



E 



{L{x + X{a^,a2),D) - L{y + X(ai, 02), D)Y 
< c" (nl)^^^'^^'^ \x — y^^"^ r^''^^~'^^"'\ 



(4.2) 



In the above, ki > is the same constant as in Lemma 



Proof For any fixed (01,02) G M , we define the Gaussian random field Y = {Y{s,t), (s,t) G 
M^} with values in by Y{s,t) = X{s,t) - ^(01,02). It follows from ([XT]) that if X has a 
local time L{x,S) on any Borel set S, then Y also has a local time L(x,S) on S and, moreover, 
L{x + X (ai , 02) , S) = L{x,S). Since X has stationary increments, both Lemma 13.41 and Lemma 
ESI hold for the Gaussian field Y. This proves (gH) and (g^]). □ 

The following lemma is a consequence of Lemma l4.ll and Chebyshev's inequality. 

Lemma 4.2 Assume the conditions of Theorem \3.S\ hold. For any 6 > 0, there exist positive and 
finite constants e G (0, 1/e), C4 3, 04 ^, €4 5, (depending on ai, 02, -^ij -^2 <i'nd d only), such that for 
all (01,02) G M^, D = OAri,Af2 (('^1) ^2)5 ''') with r G (0, e), x G M*^ and u > 1 large enough, the 
following inequalities hold: 
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(i). If Ni 7^ aid, then 



»|L(x + X(ai,a2), D) > c^.^r'^^u''-^ < exp{-bu). 



(a). If Ni = aid, then 

p|l(x + X(ai, 02), D) > c, , ^ (^g + 
(^iiij. For x, y G M"^ with \x — y\ < 1 and 7 > small, 



||L(x + X(ai,a2), D) - L{y + X{ai,a2), D) 



> C4 5 |x — yy r 



7 „/3t-k7 ,,'?t+k7 



< 



exp (— 6n). 



(4.3) 



<exp(-6n). (4.4) 



(4.5) 



Proof The proofs of Parts (i) and (iii) based on Lemma 14.11 and Chebyshev's inequality are 
standard, hence omitted. In the following we prove (ii). Define the random variable A = L(x + 
X(ai, 02), D)/r^2, For n > large, let n = \u\, the largest positive integer no bigger than u. We 
apply Chebyshev's inequality and Lemma l4.1l to obtain 



A > c u log I e + 



OL'2 OL-^ ' 



< 



ec 



r°'2-°'i 



log" e+ 



CK2 Oi-^ > 



where c > is a constant whose value will be determined later, and where we have used the fact 
that for j E {1, 2, . . . , n}, 



CV2 ^1 

jv7~ ivT 



log e H < log e H 

By taking c = c^^ large so that log {c^^/iec)^ < —b, we obtain ()4.4p . 



□ 



The following result about the exponential integrability of L{x,D) is a direct consequence of 
Lemma 14.21 We omit its proof. 

Theorem 4.3 Assume that the conditions of Theorem \3.3\ hold and let Di := OjVj^7V2(0, 1). Then 
there exists a constant 5 > 0, depending on ai, a^, Ni, N2 and d only, such that the following hold: 



(i). If Ni 7^ aid, then for every x € M'^ 

< 1, 

where r]r is the constant given in |g. i^[ ). 



(4.7) 
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(ii). If Ni = aid, then 



where ipsiy) = y/ log(e + y) for all y > 0. 

Now we study the local Holder condition of the intersection local time and its connec- 

tion to fractal properties of the set of intersection times M2 of B"^ and i?"^ . 

Since M2 is the zero-set of X, namely, M2 = X^^{0), and the Gaussian random field X satisfies 
the conditions in Xiao (2009). It follows from Theorem 7.1 in Xiao (2009) that 

dimjjM2 = dimpM2 = /3r (4.9) 

with positive probability. In the above, dimjj and dimp denote Hausdorff and packing dimension, 
respectively; see Falconer (1990) for more information. In Corollarv 14.61 below, we will show that 
(fOj) holds with probability 1. 

It is an interesting problem to determine the exact Hausdorff and packing measure functions 
for M2. For this purpose, the limsup and liminf type laws of the iterated logarithm need to be 
established, respectively for the intersection local time L°'^''^^{-). 

In the following, we consider the limsup laws of the iterated logarithm for the local time L{x,-) 
of X. By applying Lemma 14.21 [with (01,02) = (0,0)] and the Borel-Cantelli lemma, one can easily 
derive the following result: There exists a positive constant g such that for every x G R'^ and 

L{x,ON,,N,i{s,t),r)) 
mil bup J r _ L.4 Q 5 d.b. , lu J 



where 



r-- 1 loglog(l/r))''" if Ni / aid, 



r^^(loglog(l/r)) log(e+ ('°g^°g(^ff,lV ) Ni = aid. 



(4.11) 



It is worthwhile to compare (j4.10p with the corresponding results for {N, d) fractional Brownian 
motion of index a in Xiao (1997) and the (A'^, d) fractional Brownian sheets with index (a, . . . , a) € 
(0, 1)^ in Ayache, Wu and Xiao (2008). In the former case, X is isotropic and its local time L{x, •) 
satisfies 

r L{x,ONit,r)) 

hmsup^^T — T-, —r — nT7<c.y, a.s., 4.12) 

^^0 r^-"'^(loglogl/r)°'^/^ - ' ^ ^ 

while the local time of the {N, d) fractional Brownian sheet with index (a, . . . , a) G (0, 1)^ satisfies 

r L{x,OMt,r)) 

hmsup^TT — tt; — , , ■ < c, „ , a.s. (4.13) 

,.^0 r^"'"^(loglogl/r)"'=' - ^ ^ 

Note that, the anisotropy of the fractional Brownian sheet only increases the power of the correction 
factor log log 1/r. For the Gaussian random field X defined by (|1.3|) with A^i = aid, (|4.10|) suggests 
that the asymptotic properties of the local times of X may be significantly different from those in 
()4.12p and (14.130 . In fact, when A'^i = aid and as r J, 0, we have 

r^2 loglog(l/r) logloglog(l/r) if 02 = ai, ^ > j^, 
(^i(r)~<J r^^loglog(l/r) ifa2 = ai, ^<^, 

log(l/r) log log(l/r) if 02 > ai- 
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However, in this later case, it is unclear to us whether the logarithmic correction factor in (j4.1ip 
is sharp. It would be interesting to study this problem and establish sharp laws of the iterated 
logarithm for the local times of X. For such a result for the local times of a one-parameter fractional 
Brownian motion, see Baraka and Mountford (2008). 

As a consequence of (|4.10|) we have for the intersection local time of B'^^ and -B"^ that, for 
every (s,t) G M^, 

r L'-^^'''{ON,,N,{{s,t),r)) ^ 

limsup ^ -— -<c.f., a.s. (4.14) 

It follows from Fubini's theorem that, with probability one, (j4.14p holds for Ajv-almost all (s,t) G 
M^. Now we prove a stronger version of this result, which is useful in determining the exact 
Hausdorff measure of M2. 

Theorem 4.4 Assume that + > c^- -^et r e {1, 2} he the integer defined in \3.12\) and let 
D = OjVj^7V2 (0, -R) he fixed. Let L"i'"2(-) he the intersection local time of B"'^ and B"'^ , which is a 
random measure supported on the set M2 ■ Then there exists a positive and finite constant c^ g such 
that with prohability 1, 

l.msup ^"'°<""-f «'■'>■'•" <c„ (4.15) 

holds for L°'^'"^{-)- almost all (s,t) G D, where 9^1 (r) is defined in 



Proof Again we work on the random field X defined by (|1.3p . For every integer /c > 0, we 
consider the random measure Lk{x,») on the Borel subsets C of Oat^^atj (0; -^) defined by 



Lfc(x, C) = J (27rA:)^/2 exp - tl^^^El^^'^ j^^^ 




ic 



2 . (4-16) 



exp I - ^ + H^, X{s, t) -x)] d^ dsdt 



Then, by the occupation density formula (j3.ip and the continuity of the function y 1-^ L{y, C), one 
can verify that almost surely Lk{x, C) L{x, C) as k ^ 00 for every Borel set C C 0]yi,N2{^j -^)- 
For every integer m > 1, denote fm{s,t) = L[x, 0]\fj^^N2{{s,t),2~'^)Y From the proof of Theo- 
rem [331 we can see that almost surely the functions fm{s,t) are continuous and bounded. Hence 
we have almost surely, for all integers m, n > 1, 

[ [fm{s,t)r L{x, dsdt) = \im [ [fm{s,t)r Lk{x,dsdt). (4.17) 

It follows from (j4.17p . (|4.16p and the proof of Proposition 3.1 of Pitt (1978) that, for every positive 
integer n > 1, 

IE / [fm{s,t)r L{x,dsdt) 



1.^ 



^1 / / / exp ( I /^-^gx 



n+1 

X " 



Eexp I (u^X(sj,ij)) J dudt, 
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where u = (u\ . . . , m"+^) G m("+i)'^ and t = {s, t, si, h, . . . , s„, Similar to the proof of (l3T3]) 
we have that the right hand side of (j4.18p is at most 



'^4.10 "'^ 



d/2 

detCov{Xo{s,t),Xo{si,ti), . . . ,Xo(s„,t„)) 



< 



Ojvi,jv2(0,-R) JOjVj,iV2{(s,t),2-™)" 

cl^^ (n!)''- 2-™"^-, if A^i ^ aid, 

c^^^ n! 2-'^'"^2 YYJ^^ log (e + 2(°2-"lW^ if = ^^^i^ 



(4.19) 



where 04 ^^ is a positive finite constant depending on ai, 02, Ni, N2, d and R only. 

Let p > be a constant whose value will be determined later. We consider the random set 

Dm{u:) = {{S,t) G On.M^^R) fmis,t) > pM'^'")} ■ 

Denote by the restriction of the random measure L{x,-) on Otv^^atjIO) -^)) that is, Pu)[E) = 
L{x, E n Onx,N2{^i ^)) fo'^ every Borel set E CR^. Now we take n = [log mj . Then, by applying 
(j4.19p and by Stirling's formula, we have 

IE D\ \fm(s,t)]'^ L(x,dsdt) 

E.„(»,„)< '°"-''7;,,p-4 -<rrr^ (4.20) 

provided p > is chosen large enough, say, p > c^^^-^ := c^g. This implies that 

Elf^pUDm)] < 00. 

\m=l / 

Therefore, with probability 1 for //^ almost all {s,t) G Oiyi,N2i^j we derive 

.in..np£^^^?^lM£^<e.,. (4.21) 

Finally, for any r > small enough, there exists an integer m such that 2""^ < r < 2^*"+^ and 
(j4.2ip is applicable. Since <fi{r) is increasing near r = 0, (j4.15p follows from (j4.2ip . □ 

As an application of Theorem 14. 4t we derive a lower bound for the exact Hausdorff measure of 
the set M2 of intersection times. The corresponding problem for the upper bound remains open. 



Theorem 4.5 Let = {^"^(s), s G M^i} and = {^"^(t), t G M^^} ^y^^ independent 
fractional Brownian motions with values in M'^ and Hurst indices ai and 02, respectively. Assume 
that ^ + ^ > d. Then, for every R > 0, there exists a positive constant c^ such that with 
probability 1, 

ipi-m{M2nON^,N2{0,R)) > c,^,2L{x,On^,N2{0,R)), (4-22) 
where (pi-m denotes the ipi-Hausdorff measure. 

Proof As in the proof of Theorem 4.1 in Xiao (1997), (j4.22p follows from Theorem 14.41 and the 
upper density theorem of Rogers and Taylor (1961). We omit the details. □ 
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As a corollary of Theorem 14.51 we have the following result which is stronger than (j4.9p . 

Corollary 4.6 Let B"'^ and B"^ be defined as that in Theorem \4-5\ If ^ + ^ > d, then with 
probability 1, 

( N -aid if ^> d, , , 

dim,M, = dim,M, = + ^iV, - ^f"^<d<^ + ^. ^''^^^ 

Proof It is known from Theorem 7.1 in Xiao (2009) that dimpM2 < (3^ almost surely. In order 
to prove dimjjM2 > [ir almost surely, thanks to Theorem 14.51 it is sufficient to show that with 
probability 1, the intersection local time L"i'"2(O7Vi,Af2(0> -^)) > for R large enough. We can 
actually prove a stronger result than this last statement. First note that, when x = 0, (j3.17p 
becomes an equality. Thus, one can verify that E[L°i'"2(0^^^^^(o, 1))] > 0, which implies that 
L"^^''^'^[On-i,N2{^-: 1)) > with positive probability. More precisely, there exist positive constants 5i 
and 62 such that P(L"i'"2(o^^^^^(o, 1)) > 5i) > 62. 
For any integer n > 1, define the event 

L(0,[0,2-"/°i]^i X [0,2-"/"2]Af2) > 5^2""^~+^"°'H. 

By the scaling property (j3.9p . we have ¥{An) > 82 for all n > 1. It follows from this and Fatou's 
lemma that P(limsup74n,) > 62- This implies that with positive probability 

n— >oo 

lim sup , y ^ > 61 . (4.24) 

— H -—a 

Finally, note that the Gaussian field X has stationary increments and satisfies the condition of 
Theorem 2.1 of Pitt and Tran (1979), which is a zero-one law for X at 0. Hence (|4.24p holds with 
probability 1 which, in turn, implies L"i'°2(0^^_^2(o, ii)) > for all R>0. □ 



5 HausdorfF and packing dimensions of D2 



In this section, we determine the Hausdorff and packing dimensions of the set D2 of intersection 
points of and B'^^ defined hy D2 = {x e R'^ : x = B°''{s) = B'^^it) for some (s,t) G M^}. 
Note that we can rewrite D2 as D2 = B"i(M^i) n 



Theorem 5.1 Let and B"'^ be defined as that in Theorem \4-5\ If ^ + ^ > d, then with 
probability 1, 



dimjj£'2 = dimpL'2 = < 



d 

N2 

^ + ^-d 

ai a2 



if Ni > aid and N2 > a2d, 

if Ni > aid and N2 < 02^, 

if ^1 ^ ctid and N2 > a2d, 

if ^1 ^ ctid and N2 < a2d. 



(5.1) 



In order to prove Theorem lS.l^ we will make use of the following two lemmas which are corollaries 
of the results in Monrad and Pitt (1987). 



24 



Lemma 5.2 Let i?" = {B°'{t),t € R''} he a fractional Brownian motion with values in M'^ and 
index a G (0, 1). If p > ad, then almost surely i?"(MP) = M*^. 

Lemma 5.3 Let = {B^{t),t € W} be a fractional Brownian motion with values in and 
index a € (0,1). If p < ad, then for any constants R > 1, e > and /? > such that < 
a — e < P < a, the following statement holds: With probability 1, for large enough n and for all 
balls U CW^ of radius 2""^, the inverse image {B'^)^^{U) can intersect at most 2"*^^^ cubes I^j^ of 
the form 

In,-k = e [0, R]P : {ki - 1)2-" <ti< ki2-'', i = l,2,...,p}, 
where k = {ki, . . . , kp) and 1 < < i?2" for i = 1, . . . ,p. 

Proof of Theorem 15.11 We prove (jS.ip by considering the four cases separately. 

Firstly, we assume that A'^i > aid and A''2 > a2d. It follows from Lemma 15.21 that almost surely 
S"i(M^i) = 5"2(M^2) = ]^d_ jjgj^pg = 5"i(M^i) n 5"2(]^A'2) = j^_g_^ ^j^i^j^ implies that 

dimjjZ)2 = dimpD2 = d almost surely. 

Secondly, we assume that Ni > aid and N2 < a2d. Then D2 = S"i(M^i) n S"2(]R^2) = 
S"2(M^2) a.s., which yields dimjjZ)2 = dimpiD2 = N2/a2 almost surely. The proof for the case 
-^1 < did and N2 > a2d is similar. 

Finally, we consider the case of A'^i < aid and N2 < a2d [in addition to ^ + ^ > d]. Let 5*2 
be the projection of M2 on M^^. Then B°'^{S2) = D2. Since, for every e > 0, B°''^{t) satisfies a 
uniform Holder condition of order 02 — £ on every compact interval of , we have 

1 N-i No 

dimp5"2(52) < — dimp52 < — + — -d, a.s., (5.2) 

a2 ai a2 

where the last inequality follows from (j4.23p . 

It only remains to prove dimpjD2 > ^ ~^ ^ ~ ^ almost surely. For this purpose, we denote 
£ = ^ + g - d and define an (iV, 2£i)-Gaussian random field Z = {Z{s, t), {s, t) G M^} by 

Z{s, t) = (S"i (s), (t)) , (s, t) £ R^. 

Set D2 = {{x,x) : X e D2}. Then 

Z-\D2) = M2. (5.3) 

Fix an u; € such that the conclusion of Lemma 15.31 holds. Assume that for some constant 
r] > 0, dim^D2{uj) < I — r]. [We will suppress uj from now on.] Then, for any n large enough, there 
exists a sequence of balls {Ui} in R'^ with radius < 2"", such that 

D2^[JUi and ^(diamJ/i)^-'' < 1, (5.4) 

i i 

where diamC/ denotes the diameter of U. Choose positive constants e, 71 < ai and 72 < 02 such 
that 

71 < 72 and ed(^ + l) <^. (5.5) 
V71 / 2 

Let rrii and nj be integers that satisfy 

2-(™.+i)7i < diamC/i < 2-'"'Ti and 2-("'+^)t2 < diamf/^ < 2-"^T2. (5.6) 
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By (|5.3p and Lemma 15.31 we have 

M2n[0,Rf c\Jz-\u,xUi) 

' r 1 (5-7) 

C y |the union of at most 2(™»+"»)"'^ cubes /^^j x I A. 

i 

Denote the cubes in the right-hand side of (|5.7|) by Cij. Note that, since 71 < 72, we derive from 
(j5.6p that diamCjj < 3 • 2~"' for i (or n) large enough. Combining this with (j5.4p . (j5.6p and (j5.6p . 
we derive that for all n large enough, 

^^(diamC,j)^2(^-2) ^2("^«+'^0^'i(2-"«)T2(^-i) 

i j i 

<C3,5:2-(^+t)^^(2---^/-i (5.8) 

j 

< C5,3 ^(diamf/i)^"'' < c^ ^. 

i 

It follows from (jSTT]) and (fSlH]) that 

dimjM2n[0,i?]^) <72(^-|) <a2(^-|). 

Hence we have proven that, for any rj > 0, 

p|dimjjL>2 >i-v} >P{dimjj(M2n [0,i?]^) >a2(^-|)}. (5.9) 

Letting R ] 00 and 1] [ along the rational numbers and by using ()4.23p . we obtain dimjjD2 > ^ 
almost surely. This finishes the proof. □ 
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